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ABSTRACT 

We present the results of a study of the late-type spiral galaxy NGC0959, before and after applica- 
tion of the pixel-based dust extinction eorreetion described in Tamura et al. 2009 (Paper I). Galaxy 
Evolution Explorer (GALEX) far-UV (FUV) and near-UV (NUV), ground-based Vatican Advanced 
Technology Telescope (VATT) UBVR, and S'piizer /Infrared Array Camera (IRAC) 3.6, 4.5, 5.8, and 
8.0 images are studied through pixel Color-Magnitude Diagrams (pCMDs) and pixel Color-Color 
Diagrams (pCCDs). We define groups of pixels based on their distribution in a pCCD of {B — 3.6 /xm) 
versus (FUV — ?7) colors after extinction correction. In the same pCCD, we trace their locations 
before the extinction correction was applied. This shows that selecting pixel groups is not meaningful 
when using colors uncorrected for dust. We also trace the distribution of the pixel groups on a pixel 
coordinate map of the galaxy. We find that the pixel-based (two-dimensional) extinction correction 
is crucial to reveal the spatial variations in the dominant stellar population, averaged over each res- 
olution element. Different types and mixtures of stellar populations, and galaxy structures such as 
a previously unrecognized bar, become readily discernible in the extinction-corrected pCCD and as 
coherent spatial structures in the pixel coordinate map. 

Subject headings: dust, extinction — galaxies: individual (NGC0959) — galaxies: spiral — galaxies: 
stellar content — galaxies: structure 



1. INTRODUCTION 

Studying the distribution of different stellar popula- 
tions within galaxies is crucial to understand the for- 
mation and evolutionary history of galaxies. One of 
the obstacles that prevents us from directly observing 
these stellar populations is the attenuation and red- 
dening by dust that is randomly distributed among, 
or projected along the line-of-sight toward, these stel- 
lar populations. This ha s remained a major problem 
for iiearly a century (e.g.. iTrumpleil 119301: [Mathis et al.l 



1977t IV iallefond et all 1 19821: iCaplan fc Deharvene|ll985l 
Witt et al. 1992: Roussel et a l. 2005; Driver et al. 200$. 
While the main effect of dust is to attenuate the light 
emitted by stars, dust also scatters light out of the line- 
of-sight, and may re-di rect light from nearby regions into 
our l i ne-of- sight (e.g.. IWitt et al.l 119921: IWitt fc Gordo^ 
119961 12OOO0 . The impact of extinction — the wavelength- 
dependent net effect of absorption and scattering — 
correlates wi th the spatial distribution of star-f orming 
regions (e.g., IWaller et al.llT99l IDeo et al.l[2006h . as it 
depends on bo th the geomet ry of the dust distribu- 
tion (e.g..lElmc green 1980: Walt erbos fc Kennicut^l 19881: 
ICalzctti e t aLfri994 iWitt fc Go^^doillTl996l I2000D and on 
the physica l and chemical propert i es of the dust (e.g. , 
van HoutCTlfTgel IwitreTall [1991 TWhittet et ai1[200ll 
2004( 1. On tvpical observational scales, the dust will be 



intermixed with the stellar populations as thin layers, fil- 
aments, and dense clumps. Even though some light may 
be scattered into the line-of-sight from the rear and from 
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regions near the line-of-sight, most of the extinction by 
far will be due to dust that is distributed in front of the 
stellar populations of interest (|Bvunl[r992[ ). 

Popular methods to measure dust extinction, such 
as the use o f Hydrogen recombin a tion line ratios (e.g. , 
Rudyl [1983: IScoviUe et al.l [200lt ICalzetti et aH [200l 
2007t iKennicutt et al.l 120091) . 



spec - 



20091) ■ the ultraviolet (UV) , 

tral slope (e.g.. ICalzetti et al.l [l99llKong et alll2004l) . 
or ratios o f the UV and total infrar e d fluxes (e.g., 
Buat fc Xul 119961 : ICalzetti et aD 120001: iBoissier et ahl 



2004 1200511. have various constraints and limitations 
(e.g..lXu fc Heloul[l99l [Petersen fc Gammelgaa rd 1997 



ReganI I200GI: iPrice et al.l 120021: IBoissier et al.i i2004 : 



Rieke et al.ll2004D . In lTamura et al.l (|2009l . hereafter Pa- 
per I), we introduced a new method to measure the dust 
extinction in galaxies, and demonstrated that the spa- 
tial distribution of dust extinction in a late-type spiral 
galaxy, NGC 0959, can be mapped using only the flux ra- 
tio of optical V -band and mid-infrared (mid-IR) 3.6 ^bTn 
images. In the present paper, we present the results of a 
multi-wavelength (UV-optical-mid-IR) study of the ef- 
fect of the pixel-based two-dimensional extinction cor- 
rection of Paper I on color composite images, a pixel 
Color-Color Diagram (pCCD), and a pixel coordinate 
map of NGC 0959. These results show that a detailed 
pixel-based extinction correction is essential to reveal the 
nature and distribution of stellar populations in galaxies. 

2. ASSUMPTIONS ABOUT DUST EXTINCTION 

Even though the effects of dust attenuation dominate 
over those by scattering, scattering can not be ignored 
completely. IWitt fc Gordm] (|1996L l2000( ) show that, de- 
pending on dust geometry, the effect of scattering can be 
significant. Among the modeled dust geometries, their 
CLOUDY model corresponds closest to the dust geome- 
try we assumed in Paper I: a geometry where thin layers. 
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filaments, and dense clumps of dus t are intermixed with 
the stars. IWitt fc Gordoi] (|2000l their Figure 8) find 
that this model shows the smallest effect from scatter- 
ing, with fsna.t. /f* 5. 0.1 (o r A/i < 0.1 mag arcscc"^, in 
agreement with lBvunlll992| ) for all modeled wavelengths 
(A > 0.1 fim) and metallicities, even for very clumpy dust 
distribu tions. 

While IWitt fc GordonI (Il996l. [2000l) model the effect of 
dust attenuation and scattering, IGordon et aP (l200l — 
and references therein — measured the UV-near-IR ex- 
tinction toward different regions within the Milky Way 
(MW), the Large MageUanic Cloud (LMC), and the 
Small Magellanic Cloud (SMC). Since the effects of at- 
tenuation and scattering c annot be separated, the ob- 
served extinction curves in IGordon et al.l (120031 ) are the 
total effect of dust. Similarly, the extinction measured 
in Paper I represents the total, or net, effect of both at- 
tenuation and scattering. 

One needs to adopt an extinction curve to scale the 
measured visual extinction, Ay, to the extinction Ax 
in other bandpasses. There appears to be a sys tematic 
change from a n SMC ([Gordon fc Clavto n' 19981 via an 
LMC (Missch et al.lll999D t o a MW tvpe ( Clavto n et al.l 
120001 : IValencic et al.l l2003| ) extinction curve, which — 
while controversial and unproven — is often taken to re- 
flect the difference in average metallicity and/or mor- 
phology of these galaxies. Even within a single galaxy, 
however, there is significant variation between extinc- 
tion curves derived for differe nt sightlines toward indi- 
vidual stars or regions (e.g.. IGordon fc ClavtonI [19981: 
Misselt et all Il999t iClavton et al.l l2000t IValencic et al.l 
2003HGordon et al.ll2003[ ). The rapid decrease in attain- 
able spatial resolution with increasing distance renders 
similar measurements impossible for galaxies beyond the 
Local Group. Based on its overall morphology and lumi- 
nosity, we assumed in Paper I that an LMC-like extinc- 
tion cur ve — actually the "LM C2 supershell" extinction 
curve of IGordon et al.l (|2003D — would be more appropri- 
ate for NGC 0959 than either an SMC or MW-type ex- 
tinction curve. At the ~ 10 Mpc distance of NGC 0959, 
we expect that the mixing of light from neighboring 
stellar populations on ■~250pc scales will tend toward 
a mean similar to that for the LMC2 supershell, even 
when the extinction toward individual stellar populations 
might be better characterized by an SM C or MW-type 
extinction curve. Adopting the SMC bar ([Gordon et al.l 
[200l or an average MW-extinction curve, would over- 
or under-estimate the extinction at shorter wavelengths, 
respectively. 

3. DATA-SET AND PREPARATION 

The galaxy selected for thi s pilot study, NGC 0959, is a 
late-type spiral galaxy (Sdm: [de Vaucouleurs et aT 
at a distance of Z) = 9.9 ±0. 7 Mpc (I Mould et a ' 
with an inclination of 



1991) 



120001) 

'50° ([Esipov et al.lll99lD. This 



galaxy has been observed w ith GALEX ([Martin et al.l 
[2005t iMorrissev etldl [20071) in the FUV and near-UV 
(NUV), with the Vatican Advanced Tech nology Tele- 
scope (VATT) in UBVR from the ground ([Tavlor et al.l 
[2001 . and wi th Spitzer/IKAC in the 3.6, 4.5, 5.8, and 
8.0 /im filters ([Fazio et al.ll2004D . Near-IR JHKs images 
arc also a vailable from the Two Micron All Sky Survey 
(2MASS; iSkrutskie et~aI]|2006D . but these near-IR im- 
ages are too shallow for our purpose, and therefore ex- 



cluded from the analysis. All displayed images and pixel- 
maps are rotated such that North is up and East is to 
the left. The surface brightness and colors are in units of 
mag arcsec"^ and mag, respectively, and are on the AB 
magnitude system ([Okelll974l: fOke fc Gunnlll983f ). 

To perform the pixel-based analysis, all images are 
convolved and re-sampled to a matching point spread 
function (PSF) and pixel scale. Among all images 
from the different telescopes and instruments of inter- 
est, the coarsest pixel scale and PSF are 1'.'5 pixel^^ and 
- 5'.'3 FWHM, respectively, for the GALEX NUV image. 
IDL'* function "f rebin" (with its flux conserving option 
switched on) and IRAF^ routine "gauss" (with a round 
two-dimensional Gaussian convolution kernel) are used 
to match the pixel scale and PSF at all wavelengths, 
while conserving the total amount of surface flux. We 
then select only those pixels that have a signal-to-noise 
(S/N) > 3.0 in all filters for further analysis, ensuring re- 
liable pixel surface brightness measurements and colors. 
Below, we briefly describe the data used in this analysis. 
For the details of the data preparation and the method 
of extinction correction we refer the reader to Paper I. 



3.1. GALEX Images 

The GALEX FUV and NUV images are obtained 
from the Multi-Mission Archive at the Space Telescope 
Science Institute^ (MAST). NGC 0959 wa s observed in 
the G ALEX All-sky Imaging Survey (AIS: lMartin et al.l 
20051) and the G ALEX Nearby Galaxy Survey (NGS; 
Gfl de Paz et~all [20071) . Since the galaxy is barely vis- 
ible in the AIS images, wc will use only the deeper 
(~ 1695 sec) NGS data for our analysis, focusing only on 
a 7'.5 X 7.'5 image section centered on NGC 0959. Since 
the average FWHM of stars near NGC 09 59 is larger in 
the NUV image (e.g.. iMartin et al.l [20051 ). wc convolve 
the FUV image to match the PSF in the NUV. 



3.2. Ground-based Images 

Ground-ba s ed U B VR images were obtained by 
iTavlor et all ([20050 with the direct CCD imager at 
the VATT at Mt. Graham International Observatory 
(MGIO) in Arizona. These flux calibrated images 
are available through the NASA/IPAC Extragalactic 
Database^ (NED). For a de tailed descr iption of this data- 
set, we refer the reader to iTavlor et al . (2005), and ref- 
erences therein. The effective exposure times in U, B, V, 
and R are 1200, 600, 480, and 360 sec, respectively. The 
native resolution and pixel scale arc ^1'.'3 (FWHM) and 
0'.'37 pixel"^. These images are registered, convolved and 
rcsamplcd to match the orientation, PSF, and pixel scale 
of the GALEX NUV image. 

* IDL is distributed by ITT Visual Information Solutions (Re- 
search System Inc.), Boulder, Colorado: http : //rslnc . com/ldl/ 

^ IRAF is distributed by National Optical Astronomy Observa- 
tory (NOAO), which is operated by the Association of Universities 
for Research in Astronomy (AURA), Inc., under cooperative agree- 
ment with National Science Foundation (NSF): http://iraf .net/ 

6 Galaxy Evolution Explorer, GR4 Data Release (May 1, 2009): 
http : //galex . stscl . edu/GR4/ 

^ NASA/IPAC Extragalactic Database: 

http : //nedwww. ipac . caltech. edu/ 
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3.3. Spitzcr /IRAC Images 
The Spitzer /IRAC 3.6-8.0 /im pipeline-product images 
were obtained from the Spitzer Archive^ via Leopard. 
For each filter, the mosaiced image was created with 
an effective exposure time of 26.8 sec. The native pixel 
scale is 1'.'2 pixel"^, and the effective resolution (FWHM) 
ranges from ~2'.'2 at 3.6 fim to ~2'.'3 at 8.0 /im. Like for 
the VATT images, we match the orientations, PSF, and 
pixel scale of the IRAC images to those of the GALEX 
NUV image. 

4. CONSTRUCTING AN EXTINCTION MAP OF NGC 0959 

In Paper I, we estimated the visual dust extinction 
{Ay in mag arcsec"^) measured over each l'.'5xl'.'5 ~ 
72x72pc2 pixel in NGC 0959 as follows. From his- 
tograms of the observed visual to 3.6 yitm flux ratio 
(/y//3.6Aim) in each pixel, we estimate the intrinsic 
extinction-free flux ratios (/3y,o) for two groups of pixels: 
pixels apparently dominated by the light from younger 
and pixels apparently dominated by the light from older 
stellar populations. These are separable based on their 
distribution pattern in a pCMD of the observed nv ver- 
sus ([/ — 3.6/im) color. Since the mid-IR 3.6 /xm flux is 
assumed to be minimally affected by t he dust, and hence 
usually treated as e xtinction-free (e.g.. lFazio et al.ir2004l : 
iWillner et al.l 120041 ) , Ay in each pixel can be inferred 
from the difference between the observed {fv,ohs) and es- 
timated extinction-free (i.e., fv,o = Pv,o x /3.6/im,obs) V- 
band fluxes. For further details of the method, we refer 
the reader to Paper I. The extinctions in other band- 
passes {GALEX NUV and FUV, and optical U, B, and R 

do pted L 

D pool . 

In Figure [U the extinction map of NGC 0959 thus pro- 
duced, darker grayscales correspond to higher values of 
Ay- The grayscales saturate for Ay > 0.4 mag arcsec"^ 
(indicated by the white vertical line in the color bar) 
to enhance the visibility of the lower Ay values in the 
galaxy, and so differ from that of Figure 10(a) in Pa- 
per I. The maximum extinction measured in this galaxy, 
averaged over a pixel, is Avmax — 0.8mag arcsec"'^. 

5. LIFTING THE VEIL OF DUST FROM NGC 0959 
5.1. Color Composite Images 

To visually (qualitatively) investigate the effect of 
our pixel-based extinction correction on an image of 
NGC 0959, we first construct two color composites of 
the galaxy, composed of the Spitzer /IRAC 3.6/im (red 
channel), the ground-based V (green channel), and the 
GALEX FUV (blue channel) images. Figures [2ja) and 
[21(6) show the color composites before and after extinc- 
tion correction. The image resolutions are matched to 
that of the GALEX NUV image. For easy comparison, 
both images were created using the same color stretch, 
and with the IRAC 8.0 /im contours over-plotted. 

In Figure m a), there are several regions in the galaxy 
that appear much redder than other parts of the galaxy. 
These regions include: (a) the southern half of the galaxy, 
especially along the southern "edge" of the galaxy disk 
(as defined by our S'/A',„i,i ~ 3.0 requirement in each of 

® Spitzer Science Center (SCC) Data Archives/ Analysis: 
http : //ssc . spitzer . caltech . edu/archanaly 



bands) are then scaled from A y using the ad o pted _ LMC2 
supershell extinction curve of iGordon et al.l 
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Fig. 1. — Spatial distribution of the pixel- averaged visual dust 
extinction, Ay. Extinction values map onto grayscales from 
light-gray {Av = 0) to black (j4y=0.4mag arcsec~^; indicated 
by a white vertical line in the color bar). The few pixels with 
0.4 < Ay < 0.8 mag arcsec"^ are also rendered as black. Dotted 
contours trace the Spitzer /IRAC 8.0 /jm PAH emission and show 
good agreement with the visual extinction estimates, which were 
inferred from the observed l/-band and 3.6 fim fluxes only. The 
map measures ~135"xl20", has a plate scale of l''5pixel^^, and 
has North up and East to the left. 

the filters); (b) the strong 8.0 emission region that 
is running from NW to SE of the galaxy through its 
center, appearing especially redder at its northern end 
point; and (c) some localized regions at the north-western 
and eastern edge of the galaxy disk. Since dust features 
are not resolved at GALEX resolution, we cannot tell 
whether these pixels in Figure[5](a) are red due to extinc- 
tion, or are dominated by intrinsically red stellar popula- 
tions. From the distribution of the estimated visual dust 
extinction Ay (Figure [Ij, however, we infer that these 
red pixels are indeed red due to intervening dust at these 
locations. 

After applying our extinction correction. Figure 6) 
shows that many of the strikingly red pixels in Fig- 
urelUa) indeed have become bluer. Especially regions 
(b) and (c) have become much bluer than before, while 
some of the regions (a) still seem to have relatively 
red colors compared to other parts of NGC 0959's disk. 
There are some other regions which have become much 
bluer as well — mostly located in the SW region and 
northern regions of the galaxy. These regions already had 
a relatively blue hue in Figure [2]( a), but become much 
bluer and brighter after extinction correction. They 
likely are actively star-forming regions. 

To see whether these red and blue regions in Figured] 
actually correspond to dust features and SF-regions, we 
created a third color composite image with a higher spa- 
tial resolution from the ground-based UVR images at 
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NGC 0959 
(b) 



Red: 3.6 micron 
Green: V 
Blue: FUV 




AFTER Dust Corr. 



Fig. 2.— Color composite images of NGC0959 using images from Spitzer jlKKC S-G^tm (red), VATT V (green), and GALEX FUV 
(blue) at GALEX resolution (a) before and (6) after application of our pixel-based dust-extinction correction (Paper I). The Spitzer /TR,AC 
8.0 ixra emission is over-plotted as green contours. Both color composites are created using the same color stretch. Regions corresponding to 
high Ay values in Figure [l] clearly become much bluer after extinction correction (panel b). This is especially clear in the blue star-forming 
knots which generally coincide with peaks in the Spitzer 8.0 /im emission. 



NGC 0959 



VATT: UVR 
0.37'7pix 




Fig. 3.— Color composite image of NGC 0959 using VATT 
U, V, and R images at their native ground-based resolution 
(FWHM~1'.'3) and pixel scale (0"37pixel-i) and without any 
extinction correction. Some dust features and blue star forming 
regions are readily discernible. 

their native resolution and pixel scale of FWHM~1'.'3 
(matched across UVR) and 0'.'37 pixel" ^. Figure |3] shows 
that the blue regions that become bluer and brighter 
are indeed likely SF-regions. It also shows that the 
redder pixels in region (b) described above seem to be 
caused by thick dust lanes running from the NW to SE 



in the galaxy. It is hard to see whether regions (a) and 
(c) are caused by a dust lane or not, but regions (a) 
seems to be distributed around bluer clumps in Figure |3l 
The localized regions (c) do not appear to be unrelated 
background or foreground objects. Such objects usually 
have colors sufficiently different from those of the gen- 
uine galaxy pixels, that they would appear as a distinct 
and separate branch or grouping of pixels in pCMDs and 
pCCDs. Since we found no such feature in the examined 
pCMDs and pCCDs for NGC 0959 (e.g.. Figure H, we 
treat these regions as parts of the galaxy. 

5.2. A (B- 3.6 nm) versus (FUV- U) pCCD 

To study the effect of our pixel-based extinction cor- 
rection in a more quantitative way, we examined pCMDs 
and pCCDs using various combinations of images from 
FUV through 8.0 fim for significant features or groupings 
in the pixel-distribution. A distinct grouping of pixels 
in these diagrams should indicate that those pixels are 
dominated by the same or a similar mix of stellar popu- 
lations. Among different diagrams examined (not shown 
here), a pCCD of (B-3.6/im) versus (FUV-C/) color 
after extinction correction (Figure HJa)) was selected, 
because it clearly shows distinct tracks and groupings of 
pixels. In particular, distinct red and blue sequences can 
be identified. 

The uncertainties in the colors at the 25*^ and 75*^ 
percentiles are ~0.01 and ~0.02mag in (FUV — t/), 
and ~ 0.05 and ^ 0.08 mag in {B — 3.6 /im), respectively. 
Since the gap of ^ 0.2 mag at the bluer end of (FUV — U) 
color between the two sequences in Figure Ufa) is much 
larger than the photometric uncertainties in the colors, 
this separation of pixels is not a random effect. The se- 
quences are also not the result of differences in residual 
reddening, as can be seen by comparing the directions of 
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Fig. 4.— Pixel Color-Color Diagrams (pCCDs) of (B-3.6^tm) versus (FUV - C/) color for NGC 0959. Panel (a) shows the pCCD after 
our extinction correction (see Paper I) has been applied, while panel (6) shows the pCCD using the observed pixel colors. The 25"^ and 
75**^ percentile uncertainties in color are ~0.05 and ~0.08 mag in (B — 3.6/^m), and ~0.01 and ~0.02 mag in (FUV — f/), respectively. 
Reddening vectors corresponding to Ay =0.2 mag arcsec^^ are drawn in the bottom right corners. Pixel groups were selected in panel 
(a) and retain their group assignment and color-coding in panel (6). Small insets show the pCCDs, plotted over the same color ranges as 
the main panels, but omitting the color coding that might guide the reader's eye. While the pCCD before extinction correction shows no 
prominent well-separated features in color-color space, distinct sequences and pixel groupings appear after application of our extinction 
correction (panel a). 



the sequences with that of the reddening vector drawn 
in Figure |4|a). The red and blue sequences are con- 
nected in the 1.15 < (FUV — t/) < 1.65 mag color range 
by pixels with intermediate colors. At colors redder than 
(FUV — [/) ~ 1.65 mag, the pixels again somewhat sepa- 
rate into sequences with redder and bluer (5 — 3.6 yum) 
colors. 

As can be seen in various SED 



models (e.g., IBruzual fc Chariot! [200l 

Anders fc Fritze-von AlvenslebenI 120031 : iMarastonl 



20051: iKotuUa et al.ll2009n . the (FUV - U) color is very 
sensitive to age for the youngest stellar populations. 
Since most of the FUV flux is emitted by young, 
massive OB-stars, combining a GALEX and an opti- 
cal b road-band filter pro vides strong age constraints 
(e.g.. iKavirai et al.l [2007[) . The (S-3.6^m) color 
was selected empirically, while examining different 
combinations of colors. The IRAC 3.6 ^m bandpass is 
com monly used as a st ellar mass distribution tracer 
(e.g.. iWillner et ani2004[ ). because it is associated more 
with the distribution of redder and older stars. Since 
the optical B-band is generally sensitive to younger 
stars, the {B — 3.6 /zm) color can be used to distinguish 
mixtures of stellar populations with and without signifi- 
cant recent high-mass star formation. The combination 
of (FUV-J7) and (B-3.6^m) colors in Figure IH; a), 
therefore, provides a powerful diagnostic of the recent 
star formation history averaged over a pixel. 



5.3. Definition of Pixel Groups 

Based on the distinct red and blue sequences, as well 
as transition regions identified in Figure |4| a), we sepa- 
rated pixels into six different groups (color-coded in Fig- 
ure Ufa)). The group boundaries (selection criteria) are 
summarized in Table 1. The exact location of the bound- 
aries between some of these pixel groups is somewhat ar- 
bitrary, but is motivated by the following considerations. 

The Group I pixels (color-coded purple) reside on the 
blue part of the blue sequence, and have blue colors in 
both {B - 3.6 ^m) and (FUV - U). The dominant stellar 
populations are very young and massive OB-stars, while 
older stellar populations contribute little to the total flux 
in these pixels. The (i? — 3.6^m) and (FUV — J7) colors 
become redder as these young stellar populations age and 
the number of remaining OB-stars decreases. 

On the other extreme, the Group VI pixels (color- 
coded red) have red colors in both (i? — 3.6 /im) and 
(FUV ~ U). These pixels are therefore dominated by old, 
quiescent, "red-and-dead" stellar populations and do not 
contain a detectable fraction of young, massive stars. 

The Group II pixels (color-coded blue) have the same 
blue (FUV — U) color range as Group I, but redder 
(B — 3.6/im) colors. The blue (FUV— ?7) color implies 
the presence of OB-stars. The redder {B — 3.6 fim) color 
indicates that Group II pixels have a non-negligible con- 
tribution to the total light from older (underlying or su- 
perposed along the line-of-sight) stellar populations. 

Group V pixels (color-coded orange) cover the same. 
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TABLE 1 

Group Boundaries in the {B - 3.6 fiu) versus (FUV — C/) pCCD of Figure [H a) 



Pixel Group Color Code 



Group Boundary Definition (mag) 



Group I Purple 
Group II Blue 
Group III Green 



Group IV 
Group V 
Group VI 



Yellow 



Orange 



Red 



(FUV-C/) < 1.15 
(B-3.6/im) < 0.188 {FUV-C/) - 0.118 

(FVY-U) < 1.15 
(B-3.6/im) > 0.188-(FUV-C/) - 0.118 

1.15 < (FUV-C/) < 1.65 
(B-3.6/i.m) < 0.188-(FUV-(7) - 0.118 

1.15 < (FUV-C/) < 1.65 
(B-3.6/i.m) > 0.188-(FUV-C/) - 0.118 

(FUV-C/) > 1.65 
(B-3.6Atm) < 0.22 and {B-3.6fim) < 0.188(FUV - C/) - 0.118 

(FUV-C/) > 1.65 
(B-3.6^tm) > 0.22 or (B-3.6Atm) > 0.188(FUV - C/) - 0.118 



red, (FUV — U) color range as Group VI, indicating that 
they too lack a detectable fraction of young, massive OB- 
stars. Yet, their bluer (i? — 3.6//m) color suggests that 
the flux in these pixels is dominated by light from inter- 
mediate age stellar populations. 

The Group III and Group IV pixels (color-coded green 
and yellow, respectively) are located between these ex- 
treme cases. The light in these pixels is likely dominated 
by stellar populations in transition between the extreme 
groups along either blue or red sequences, or between 
the blue and the red sequence, or represents a mixture of 
stellar populations (from unresolved adjacent regions or 
regions superposed along the line-of-sight) with different 
star formation histories. 

5.4. Effect of Pixel-Based Extinction Correction 

Figure 131^6) shows the pCCD of (B ~ 3.6 fim) versus 
(FUV — U) color before the application of the extinction 
correction. For each pixel, the color coding was pre- 
served from that in Figure |4l[ a). This allows us to track 
the effect of applying the extinction correction in this 
color-color space. Interestingly, as shown in the inset of 
Figure mj?)), once the color coding is removed, obvious 
features like the blue and red sequences of Figure |4l( a) 
are no longer discernible. The Group I and II pixels, for 
example, which were clearly separated in Figure |4l[a), 
have blended to form a continuous distribution in Fig- 
ure \M,b). The main difference between Groups I and 
II, which both arc characterized by the presence of OB- 
starS; is the fraction of light contributed by older stellar 
populations. Group III and IV pixels, defined as hav- 
ing 1.15 < (FUV — ?7)< 1.65 mag after extinction cor- 
rection, are affected differently by the extinction correc- 
tion (as a comparison of Figures Sl^a) and|3l^fe) shows). 
Like Group I and II pixels. Group IV pixels are found 
scattered out to much redder (S — 3.6/im) colors before 
extinction correction. Most of the Group III pixels, on 
the other hand, can be located in the same color-color 
space in both panels of Figure 2] The same difference 
is seen between pixel groups V and VI. This indicates 
that Group III and V pixels have almost no measurable 
dust extinction, while Group IV and VI pixels are signif- 



icantly obscured by dust. As a side note, the reduction 
of the scatter going from the observed to the extinction- 
corrected version of the {B — 3.6 /im) versus (FUV — U) 
pCCD is only possible if the extinction correction applied 
to each individual pixel was appropriate, at least to first 
order. 

The pixel-averaged dust extinction estimated by 
our method in NGC 0959 (Figure [T]) indeed varies 
from Av~Q.O to Ay ~ 0.8 mag arcsec"^. In Pa- 
per I, we also measured the average visual dust ex- 



+0.086 
-0.049 



tinction across the entire galaxy. Ay = 0.064 
mag arcsec"^. and the azimuthally averaged radial 
extinction profile, Ay(R), which has a central value 
of Ay(R = 0)~0.25 mag arcsec"^, but quickly drops 
to below Ay = 0.1 mag arcsec"^ beyond a radius 
of '--^0.2R25, where R25 denotes the major axis ra- 
dius at the ur = 25.0 ma g arcsec"^ isophote (RC3; 
Ide Vaucouleurs et al.lll99il ). These average Ay values 
are much smaller than the peak pixel-based Ay estimate. 
Also, while the Ay and radial profile do not become 
zero at any radius, with our method (see Paper I) about 
~ 55% of the analyzed pixels in NGC 0959 have Ay = 0.0 
mag arcsec"^ to within the photometric uncertainties. 

Before we proceed to the next section, let us con- 
sider what might cause the difference in dust extinc- 
tion between Groups I and III on the blue sequence. 
Dust and gas are easily reinoved by sta rburst heat- 
ing (e.g., iMihos fc Hernguistl 1 1994 119960 and stellar 
winds (e.g. iMurrav et al' 1 120051) . and both mechanisms 
become stronger as the size of an OB association be- 
comes larger. Therefore, if the dominant stellar popu- 
lations in Group III pixels are indeed evolved Group I 
populations wherein OB-stars have died off, we would 
expect them to suffer less extinction, unless some mech- 
anism for rapid reformation and pervasive distribution 
of dust were at play. Figures |4l( a) and|H6) show that, 
once dust-free, these populations do not become signifi- 
cantly dustier evolving from Group III to Group V, ap- 
parently ruling out such rapid reformation. Since we do 
not observe a similar strong drop in extinction going from 
Group II to Group IV, we conclude that most of the ex- 
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tinction in those pixels is not physically associated with 
the OB-associations themselves, although they could still 
represent star formation "blisters" on the far side of large 
molecular cloud complexes. 

5.5. Is NGC0959 Unique? Application to Other 
Galaxies 

To determine whether our pixel-based extinction cor- 
rection can be used as a general method, or whether 
NGC 0959 is a special case, we applied the same 
method to NGC 7320 (SA(s)d at = 14.0 ± 1.0 Mpc, 
sampling 102xl02pc2 pixel"!) .^^^ UGC 10445 (SBc 
at 1? = 20.0 ± 1.4 Mpc, 146 x 146 pc^ pixel"!), which, al- 
though more distant, are selected from our larger sample 
of 45 galaxies with FUV through mid-IR imagery as rel- 
atively close analogs. A more detailed analysis of the 
distribution of extinction and the intrinsic stellar pop- 
ulations of these galaxies will be performed in Tamura 
et al. (2010; Paper III, in preparation). Here we sim- 
ply compare their (i? — 3.6//m) versus (FUV — C/) pC- 
CDs before and after the application of our pixel-based 
extinction correction, to help validate the current results 
for NGC 0959. 

Figure [5] shows the distribution of the pixels of 
NGC 7320 and UGC 10445 that meet our S/N criteria in 
the same color-color space as Figured The main panels 
show the distributions after, the insets show the distri- 
butions before extinction correction. For each galaxy, 
both the main panel and inset show the same ranges in 
color. The vertical dotted lines in the insets represent 
(FUV — [/) = 0.0 mag. Reddening vectors corresponding 
to Ay = 0.2 mag arcsec"^ are drawn in the lower right 
corners of the main panels. Since the area subtended by a 
single pixel is larger than was the case for NGC 0959, the 
blending of light from distinct stellar populations within 
a single pixel becomes more significant, resulting in less 
clear separation of blue and red sequences. Yet, different 
groupings of pixels, analogous to those defined in Fig- 
ure|ll[a) for NGC 0959, are still recognizable in Figure [5] 
for both galaxies. 

5.6. A Pixel Coordinate Map of Stellar Populations 
within NGC 0959 

Having shown that application of our pixel-based 
extinction correction reveals significant groupings of 
pixels — dominated by different types or mixtures of stel- 
lar populations — in the (S — 3.6/im) versus (FUV — J7) 
pCCD, our next question becomes whether the spatial 
distribution of pixels belonging to these pixel groups will 
reveal meaningful large- and small-scale physical struc- 
tures within NGC 0959. To address this question, and to 
study how each pixel group relates to the visual and phys- 
ical properties of stellar populations, we plot the pixel 
groups defined in Figure ^a) onto a two-dimensional 
pixel coordinate map (Figure IH]). Each square in Fig- 
ure [6] represents a 1'.'5 X 1'.'5 pixel, color-coded according 
to the pixel group it is a member of. Black dotted con- 
tours again trace the Spitzer /IRAG 8.0 /im PAH emis- 
sion. We emphasize that each pixel group was defined 
without using any pixel coordinate information. In other 
words, this is the first time that the spatial distribution 
of pixels belonging to each of the selected pixel groups is 
revealed. 



First and foremost, Figure [S] demonstrates that the dif- 
ferent pixel groups are not distributed randomly across 
the face of NGC 0959's galactic disk, but cluster in well- 
defined regions. No systematics in the data process- 
ing and process of defining the pixel groups would be 
expected to result in spatial artifacts larger than 3x3 
pixels. Since most regions cover much larger contigu- 
ous areas, we conclude that also the smallest spatial 
groupings must be genuine. Since defining the differ- 
ent pixel groups was possible only in the extinction cor- 
rected (S-3.6/im) versus (FUV - C/) pCCD, a two- 
dimensional extinction correction is thus crucial to re- 
veal regions with systematically different star-formation 
histories. 

Group I pixels (purple) are found in and around some 
of the bluest regions in the extinction-corrected color 
composite of FUV, V, and 3.6 /xm images (Figure [21(6)), 
most of which are recognizable as such in the uncor- 
rected, higher resolution, ground-based UVR color com- 
posite (Figure [3]) as weU. Since the stellar populations 
dominating these pixels are much brighter than those in 
other regions (Figure [3]), these regions must indeed be 
vigorously star-forming. 

Group II pixels (blue) are generally distributed in re- 
gions adjacent to the Group I pixels. Each of the re- 
gions (b) and (c) of Section 15.11 are found to belong 
to this pixel group, as well. In particular, much of a 
more or less linear structure that runs from the NW to 
SE of the galaxy through its center, with high 8.0 /im 
surface brightness (PAH emission) and signs of higher 
than average obscuration (Figures [1] and [3]) , consists of 
pixels that belong to Group II. These regions have the 
blue (FUV — U) colors characteristic of OB associations, 
but are redder in [B — 3.6 //m), because a larger fraction 
of the light is contributed by older stellar populations. 
These observations are consistent with the presence of 
a stellar bar. Although NGC 0959 has been stu died be- 
fore (e.g.. lEsipov et al.lll991l: iTavlor et alll2005[ ). to our 
knowledge this bar has not previously been reported. 
The dust-corrected color composite in Figure [2^6) also 
suggests that NGC 0959 has a non-negligible bulge or 
central condensation, that is partially obscured by dust 
in Figures ^a) and [31 We suggest, therefore, that its 
morphological classification be changed fro m Sdm (RC3 ; 
Ide Vaucouleurs et al.ll99il) or Sc/Irr fUGC: lNiIsMJIl973l) 
to SBcd. 

An interesting contrast in average stellar population 
age is revealed by comparing the spatial distributions 
of Group HI (green). Group IV (yellow) and Group V 
and VI (orange and red) pixels. Whereas Group HI 
pixels arc mostly found in the northwestern half of the 
galaxy. Group V and VI pixels are distributed predomi- 
nantly toward the eastern and southeastern periphery of 
the galaxy (although a smaller region of Group V pixels 
appears along the northern rim). The regions occupied 
by Group III pixels appear neither particularly blue (ac- 
tively star- forming) nor red (quiescent) in Figure [2^6), 
consistent with the idea that intermediate age (a few 100 
Myr) stellar populations are the dominant contributors 
to the flux in these pixels. The extinction map (Figure[T]) 
shows that these pixels suffer no, or at most minimal ex- 
tinction by dust. The Group V regions correspond to 
regions that have a smooth appearance without abrupt 
changes in surface brightness in Figures [2] and [3l and 
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Fig. 5.— Extinction corrected pCCDs for NGC 7320 {SA(s)d, Z) = 14.0 ± 1.0 Mpc) and UGC 10445 (SBc, D = 20.0 ± 1.4 Mpc), two more 
distant analogs of NGC 0959 (D = 9.9 ± 0.7 Mpc). Since the area subtended by a single pixel is (much) larger in these galaxies, the effects 
of light blending are more severe, resulting in smaller and fuzzier separations of red and blue sequences than for NGC 0959. Of course, 
their star formation histories may also differ. Yet, groupings of pixels that are qualitatively similar to those in NGC 0959 are recognizable 
in both pCCDs. The insets, plotted over the same color ranges as the main panels for each galaxy, show the distribution of pixels before 
dust correction. The vertical dotted lines in the inset panels indicate the locations of (FUV — U) = 0.0 mag, for comparison with the main 
panels. 



lacking any signatures of dust (Figures [T] and [3]) . Their 
red (FUV — U) colors and neutral {B — 3.6 ^m) colors in- 
dicate that these must be regions with intermediate age 
(^1-2 Gyr), largely unattenuated stellar populations in 
the outskirts of the galaxy. The Group VI pixels (red) in 
the southern half of the galaxy disk suffer some extinction 
(Figurc[T]), but the low surface brightness in these regions 
does not allow one to easily discern dust features in Fig- 
ures [2] and [3l Even though these regions suffer some ex- 
tinction, Figure lU^ a) demonstrates that these pixels are 
red mostly because the light is dominated by older (likely 
older than a few Gyr) stellar populations. The Group IV 
pixels are distributed throughout the galaxy between the 
other populations, but perhaps mostly surrounding the 
Group II pixels. Like Group III, their (FUV — U) col- 
ors indicate the presence of relatively young stars and 
absence of OB-stars, but the redder {B — 3.6 fim) colors 
of Group IV pixels show that a larger fraction of the 
flux is contributed by older populations. The extinction 
map (Figure [T]) shows that these pixels tend to suffer 
somewhat less attenuation than Group I and II pixels. 
The overall impression is that of a large-scale past star- 
formation episode that started in the south or southeast 
and propagated toward the northwest, and that is un- 
related to the ongoing massive star formation traced by 
Group I and II pixels, which presently is concentrated 
more toward the western half of NGC 0959. 

The apparent ridge of somewhat higher extinction 
along the southeastern edge (as defined by our S/N > 3 
criterion) of NGC 0959, so striking in the extinction map 



(Figure [T|) , might hint at the presence of either a warp 
in the galactic disk or an outer spiral arm delineated by 
a dust lane. In either case, no significant star-formation 
must have been associated with that portion of the galac- 
tic disk for the past few Gyr. 

Concerning the SE-to-NW stellar population gradient 
of NGC 0959 and its SE extinction ridge, it is worth 
noting that NGC 0959 is a member of the NGC 1023 
galaxy grou p, which contains 12 other major galaxies 
(|l\IIIvlll98l . The center of that group is located ~1.6° 
(-240 kpc at D~10Mpc) south of NGC 0959. The 
group environment and/or gravitational interaction with 
NGC 0949 — the group member nearest in projection on 
the sky — may well have played a role in producing the 
extinction ridge and the stellar population gradient in 
the disk of NGC 0959. We believe this certainly deserves 
further study. 

6. SUMMARY AND CONCLUSIONS 

We have presented the results of a study of NGC 0959 
using color composite images, a pCCD, and a pixel coor- 
dinate map, and demonstrated the importance and po- 
tential of a pixel-based extinction correction. Our study 
combined ground- and space-based surface photome- 
try, ranging in wavelength from the far-UV (GALEX) 
through niid-IR {Spitzer /IRAC). Among the possible 
combinations of color-magnitude and color-color dia- 
grams examined, we found that the (B — 3.6/xm) versus 
(FUV — U) pCCD proved the most powerful diagnostic 
of differences in the stellar populations contributing to 
the flux in a pixel. After applying the pixel-based, two- 
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Fig. 6. — Pixel coordinate map of NGC 0959, showing the spa- 
tial distribution of pixels belonging to the different pixel groups 
selected in Figure [Ij a), i.e., after application of our pixel-based 
extinction correction. Each square represents a single l'.'5xl''5 
pixel. The black dotted contours trace the Spitzer /IRAC 8.0 
emission. Group I (purple) and Group II (blue) pixels appear to 
define some large-scale structure in the galactic disk, highlighted 
by OB-associations. Especially, Group II pixels reveal a previously 
unrecognized bar-like structure, running from NW to SE across the 
galaxy center. 

dimensional extinction correction described in Paper I, 
we defined six different pixel groups in this diagnostic 
diagram, ranging from pixels for which both (FUV — U) 
and {B — 3.6 iim) colors indicate that their fluxes are 
dominated by very young, massive stellar populations, 
to pixels that appear to sample only light from evolved 
stellar populations. We demonstrated that it was not 
possible to meaningfully define such pixel groups before 
extinction correction. We then showed that pixels that 
belong to a given pixel group form well-defined, contigu- 
ous regions in a pixel coordinate map, revealing system- 
atic spatial variations in the dominant stellar popula- 
tions that would not be readily discernible without a two- 
dimensional correction for extinction. We were able to 
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report the presence of a previously unrecognized stellar 
bar, for example. 

Our pixel-based two-dimensional method to correct for 
extinction, based on only 3.6 fiui and y-band images 
(adding U improves fidelity), has the advantage that it 
is applicable to any galaxy that is significantly resolved 
at rest-frame ^^3.6 fim and that it is not restricted to a 
few individual sightlines within a galaxy, nor to the very 
nearest galaxies. Although seemingly simple and crude, 
we have demonstrated that application of this method 
allowed us to uncover relatively detailed spatial informa- 
tion on the nature of stellar populations and on the star 
formation history within NGC 0959, despite significant 
spatial variations in the attenuation by intervening dust. 

While the spatial resolution in the present study was 
hmited by the GALEX NUV PSF of 5'.'3 (FWHM) to 
linear scales of ^250 pc, in the near future, combination 
of HST /WFC3 and JWST will provide much higher res- 
olution (FWHM~0'.'06-0'.'ll) at rest-frame UV through 
mid-IR wavelengths. These instruments will thus al- 
low us to study the distribution of dust extinction and 
the underlying stellar populations at high fidelity and 
at much finer spatial resolution within nearby galaxies, 
as well as at similar scales of a few hundred pc within 
more distant galaxies. Since galaxies at z > 1 are impor- 
tant building blocks of the galaxies that we see today at 
z ~ 0, studying these higher redshift objects with proper 
two-dimensional extinction estimates will be of signifi- 
cant importance to reveal the nature of galaxy assembly 
and evolution. 
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